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Abstract. We study orthogonality relations for Fourier frequencies and complex exponentials in Hilbert 
spaces L 2 ([i) with measures [i arising from iterated function systems (IFS). This includes equilibrium mea- 
sures in complex dynamics. Motivated by applications, we draw parallels between analysis of fractal measures 
on the one hand, and the geometry of wavelets on the other. We are motivated by spectral theory for com- 
muting partial differential operators and related duality notions. While stated initially for bounded and 
open regions in M. d , they have since found reformulations in the theory of fractals and wavelets. We include 
a historical sketch with questions from early operator theory. 
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1. Introduction 

In this paper we give an overview (supplemented with new results and open problems) of recent research 
on orthogonality relations for Fourier frequencies, or rather the study of complex exponentials, see details 
below. This will be used in drawing up parallels between analysis of fractal measures on the one hand, and 
the geometry of wavelets on the other. 

This allows us to further include new results here involving non-linear complex systems. 

Much of our motivation comes from a problem dating back to the use of von Neumann's Spectral The- 
orem in the study of commuting partial differential operators. This led to a formulation of a number of 
duality notions. While they started for bounded and open regions in R d , many of the questions have found 
reformulations in the theory of fractals and wavelets. At the end of our paper, we include a historical section 
making connections to questions from the beginning of operator theory. 
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It is the purpose of our paper to draw attention to these interconnections. For the benefit of the reader, 
we composed in the Appendix a personal and historical section dealing with spectrum-tile duality the way 
it was first formulated by I.E. Segal and B. Fuglede. 

Specifically, we study Z/ 2 (/^)-Hilbert spaces with an emphasis on measures /i arising as equilibrium measures 
for iterated function systems (IFSs) in R d , i.e., for finite contractive function systems in R d , so in d real 
dimensions. 

One of the applications of these IFSs, and their spectral theory, is to image processing, see e.g., |Bar06j 
and |DJ06aj . The viewpoint in Bar06] is that a variety of images, e.g., large bit-size satellite photographs, 
may be condensed into IFS-codes. The trick is that such large-scale images can be shown to allow effective 
IFS-representations; and that the IFSs are encoded with a very small set of parameters. This of course is a 
feature that IFSs share with wavelets, see also [DJ06a . 

The paper is structured as follows: in Section [5] we present the parallel between multiresolution theory 
and the theory of orthonormal Fourier series on fractals. In Section [3] we formulate several open problems 
related mainly to the fractal context. In the Appendix we describe the history and motivation of the subject. 
At several points in the paper we will make connections to iteration of complex rational maps. 

1.1. Complex systems. While affine iterated function systems (IFSs) are traditionally defined for a fixed 
and finite set of transformations S = (t.;) in R d , it is often helpful to realize them by embeddings into 
compact quotients. Specifically, starting with a particular system 5, to get a fractal iteration limit X (see 
Definitions 12.21 and I2.9|) , one usually imposes suitable contractivity properties on the maps Ti making up a 
particular system S. If the mappings are affine it is occasionally possible to embed an associated fractal 
X in a compact torus T d = M. d /T where T is a suitably chosen lattice in K d , e.g., T = Z d . 

This makes an intriguing connection to complex dynamics | Bea 91, Bea92 | and Hilbcrt spaces of analytic 
functions which we address below, and again in section 3. But the connection to compact subsets in the 
complex plane opens the door to the use of Hubert space bases in the study of Julia sets derived from the 
iteration of polynomial mappings z >— » R(z) in the complex plane, see }Bea91j . 

A main focus here will be on affine fractal measures [A, but also including measures fi which are supported 
on Julia sets for certain rational mappings of one complex variable. 

In the case of measures fi supported on compact subsets X in the complex plane C, orthogonality relations 
for Fourier frequencies z n may be studied in terms of the moments, i.e., M n (/z) := J z n d\i(z) for n = 0, 1, . . . . 

The following example illustrates both orthogonality and linear relations within the set of monomials 
z n in L 2 (p). Our discussion includes measures /x; both when L 2 (p) is finite and infinite-dimensional. The 
infinite-dimensional cases often arise as IFS-limits. 

Indeed, this yields attractive classes of measures ju with compact support in the complex plane C such that 
L 2 ((i) is infinite-dimensional, and such that there are non-trivial linear relations governing the monomials 
{z n | n € No} restricted to the support of /i: Specifically, when each z n is viewed as a vector in the Hilbert 
space L 2 (p). 

This may be seen by a choice of certain fractal measures /x having compact support in C, and chosen in 
such a manner that a family of Gram-Schmidt polynomials Pk{z), k = 0, 1, 2, . . . are in fact monomials z n , 
but with n restricted to an infinite lacunary progression in No; a rather sparse subset of No- 

Using a construction of Jorgensen-Pedersen |JP98j it is possible to get the lacunary progression z n making 
up an ONB in L 2 (/i) for values of n arranged to be sparsely distributed in No in such a manner that gaps-sizes 
are powers of 4. 

We refer to |JP98j for a complete discussion of this class of examples. But for the benefit of readers, a 
summary is included here. The details of the construction further helps make the connection between the 
setting of |JP98j and the present somewhat different set of questions. 

Details: Arrange a 27r-rescaling of the circumference of the circle T in C, so it scales to one. Now construct 
an IFS measure /x with a certain compact and fractal support X contained in T as follows: Subdivide the 
circumference into 4 equal parts, and define the two mappings {tq,ti} that select every second of the 4 
sub-intervals. In a linear scale, we may take Tq(x) — x/A, and t\(x) — (x + 2)/4. Pick probabilities (1/2, 1/2). 
Then by Hutchinson's theorem |Hut81j . this system has a unique equilibrium measure fi, and the theorem 
in |JP98j states that the complex monomials {z n \ n = 0,1,4,5, 16, 17, 20, 21, . . . } form an ONB in L 2 (p). 

The following curious conclusion follows: 
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Corollary 1.1. Every function f in L 2 (n), for \x the scale-4 measure, is the restriction to the scale-4 Cantor 
set X — supp(fi) of a Hardy H 2 -function, say F . 

Specifically, we identify F with an analytic function on {z \ \z\ < 1}, and with its a.e. non-tangential limit 
to the boundary circle T (the one-torus) . And we further identify X with a subset of T as described above. 

Remark 1.2. The corresponding result is clearly not true for L 2 (T,Haar), i.e., if the Cantor measure \i 
on T is replaced with the Haar measure; or equivalently with the restricted Lebesgue measure to a period 
interval. 

In contrast, it was proved in |JP98j that the analogous construction with the middle-third Cantor set 
yields no L 2 (^t)-orthogonal subsets of {z n \ n £ No} of cardinality more than two; in fact there cannot be 
more than two orthogonal Fourier frequencies, even allowing real (non-integer) frequencies. 

In the scale-4 fractal case, it can be shown |JP98| that the Hausdorff dimension of \i and its support X 
is 1/2. Hence, as a subset of T, X is "thin", but still with the cardinality of the continuum. 

Orthogonality: The sparse distribution in No of the exponents in the ONB for L 2 (fx) for the scale-4 case 
is as follows: The values of n which turn the monomials z n into an ONB in L 2 (n) are represented by the 
following set A of all finite sums A = jo + j\4 + + J34 3 + . . . where each is a — 1 bit, i.e., each j). is 
chosen in {0, 1}. 

For related constructions, see also [D J06a|, \D J07a[ IT) J07b] . 

Corollarv ll.il is a definite statement about the specific polynomial zhz 4 and its iterations. It relates 
to functions F(z) that are analytic in the interior of the filled Julia set; in this case the closed disk {z £ 



Definition 1.3. Let R = R(z) be a polynomial. We will view R as an endomorphism in the complex plane, 
or in the one-point compactification, the Riemann sphere. The n-fold iteration of z 1— > R(z) will be denoted 
R n , i.e., inductively, R n+1 {z) := R(R n {z)), for n = 1, 2, . . . . The filled Julia set Y(R) of R consists of the 
points with bounded orbits, i.e., Y(R) :— {z \ R n (z) is bounded }. The Fatou set F(R) is the largest open 
set in which the sequence (R n ) is normal, and the Julia set J(R) is the complement of F(R). It is known 
Bea9j] that J(R) is the boundary of the filled Julia set Y(R), when Y(R) is not empty. 

It would be interesting to prove an analogous result for z 1— > z 2 + c, or for z > R(z) where R is some 
polynomial. Then the question below will be about functions F(z) which are analytic in the interior of the 
filled Julia set for the particular polynomial z 1— » R(z). 

Recall, in general, the boundary of the filled Julia set for R is the Julia set J(= J(R)) coming from the 
dynamics of R; and moreover that the IFS-Brolin measure \xr is supported on J(R). 

Let z 1 — ► R(z) be given, and let's introduce the Brolin measure [Bro65] hr on J(R). Let H 2 {J{R)) be the 
functions in L 2 (J(R), [ir) which have analytic continuation to the interior of the filled Julia set of R. 

Definition 1.4. Let z 1— > R(z) be a given polynomial. Let P be a finite set of distinct polynomials each of 
degree less than the degree of R. Let K 2 (R, P) be a space of functions F in H 2 (J{R)). We say that F is in 
K 2 (R, P) iff there are functions F p (z) in H 2 (J(R)) indexed by P such that 



It is easy to see that each K 2 (R,P) is a closed subspace in H 2 {J{R)), so in particular it is a Hilbert space. 

Remark 1.5. The finite family P that enters (jl.ip is a family of generalized filters. It depends on the 
particular polynomial z t— > R(z), and we can expect to find solutions to (jl.ip in Definition 11.41 from the kind 
of representations of the Cuntz algebras the authors studied in their work on wavelets on fractals, |D J06a| . 

A refinement of the conclusion in Corollary 11.11 is this: For R(z) := z 4 , set P := {1,2}. Hence the 
functions F in K 2 (R,P) have the representation F(z) — Fo(z 4 ) + zFi(z A ) for a pair of functions Fq and 
F% in the Hardy space H 2 . Our Corollarv ll.il above then states that there is an isometric isomorphism of 
L 2 (J(R),fi R ) with the Hilbert space K 2 (R,P). 

To illustrate Corollary 1 1.1 1 and (|1.1| . we supply the proof details in the case of R(z) := z 4 , and P — {1, z}. 



C||*|<1}. 
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Since the spectrum for \i is the set A = {^£_ Jfc4* \ ,jk & {0, 1}} C No, the family {e^ | A € A} is an 
ONB in L 2 (/.i) = L 2 (X, fi) with X C T, and we get the following L 2 (/it)-orthogonal expansion / n Ff 

(1.2) F f (z) = J2^ x 

AeA 

valid for all / G L 2 (/j,). By this we mean that every / G L 2 (fi) is the boundary representation of the function 
F f (z) in (fT5|) . analytic in the disk L> := {z G C | \z\ < 1}. 

Specifically, (|1.2|) is absolutely convergent in Z); and applying Schwarz to p. 21) we get 

l^/WI < ||/|U» M -— L=, (zeD). 
v 1 - M 

To get the desired representation 

(1.3) F f (z)=F (z 4 )+zF 1 (z 4 ) 

from Definition EH note first that A = 4A U (1 + 4A). 
Substituting this into (|1.2p then yields 

F f (z) = J2 c iXl (^) Xl +^J2 c 4A 2+ i(zY 2 , 

AiGA A 2 GA 

which is the desired formula (|1.3[) . 

Since the functions e\(x) — z x form an ONB in L 2 (fi), the representation in (|1.3[) is an orthogonal 
expansion, i.e., we have 

H/llLo) = ll^/llffa = ll^o||ff2 + H-Flllfl-2- 
2. MULTIRESOLUTION WAVELETS AND ITERATED FUNCTION SYSTEMS 

The purpose of this section is to present some key results regarding geometry and spectral theory for 
wavelets and for the fractal measures which arise from iterated function systems (IFSs). Our focus is on 
those wavelets which come from multiresolutions (MRAs), see subsection 12.31 below. This includes wavelets 
on fractals, as presented in the paper |DJ06aj by the co-authors. 

As mentioned in the Introduction, the IFSs include dynamical systems defined from a finite set of affine 
and contractive mappings in R d , or from the branches of inverses of complex polynomials, or of rational 
mappings in the complex plane. 

In terms of signal processing, what the two have in common, wavelets and IFSs, is that large scale data 
may be compressed into a few functions or parameters. In the case of IFSs, only a few matrix entries are 
needed, and a finite set of vectors in R d must be prescribed. As is shown in [Bar06], this can be turned 
into effective codes for large images. Similarly (see |Jor06j ) discrete wavelet algorithms can be applied to 
digital images and to data mining. The efficiency in these application lies in the same fact: The wavelets 
may be represented and determined by a small set of parameters; a choice of scaling matrix and of masking 
coefficients, i.e., the coefficients (a/.) in the scaling identity (|2.5p below. 

The main characters in the chapter are attractors (section l2.ip understood in the general sense so they 
include both wavelets and IFSs. In section |2"?21 we introduce multiresolution analyses (MRAs) in the form 
of scales of nested subspaces in a fixed Hilbert space, which may be L 2 (R d ), or some other Hilbert space 
derived from a computational setup, for example for iteration of rational functions in the complex plane 
(subsection II. ip . 

Our scaling identity is discussed in section 12. 3| but again in a setting which includes both the traditional 
wavelets, and their fractal variants. In subsections 12.41 - 12. 6[ we study orthogonality relations from the view 
point of a transfer operator; an infinite-dimensional operator analogous to the Perron- Frobenius operator 
from matrix analysis . 

It turns out that there is a Perron-Frobenius operator in wavelet and fractal theory which encodes or- 
thogonality relations. We call it "the transfer operator", or "the Perron- Frobenius-Ruelle operator", but 
other names may reasonably be associated with it, see e.g., |BJ02| . In the context of wavelets, the operator 
was studied in |Law91j . It will be specified by a chosen filter function too, and we denote the associated 
operator by R mo . 

Sections 12.71 and 12.81 are about applications of the transfer operator. The term "super- wavelet" in section 
12.81 refers to a super-structure. It is proved that non-orthogonal wavelets, satisfying a degree of stability 



FOURIER SERIES ON FRACTALS 



5 



(Parseval frames), can be turned into orthogonal wavelets in an enlarged Hilbert space. This is done in 
section 12.81 Motivations are from signal/image processing. By removing correlations between frequency 
bands, or between nearby pixels, redundancies are removed, and reduction in data size, and data storage, 
can thereby be achieved. 

Definition 2.1. For A G R d , let e A (x) = e ^ lXx ', x G R d . A measure fj, on M. d is said to be a spectral measure 
if there exists a set A C M. d such that the family {e\ | A G A} forms an orthonormal basis for L 2 (fi). In this 
case, the set A is called a spectrum of the measure /x. 

We are interested in the question: what spectral measures can be constructed, other than the Lebesgue 
measure on an interval? A surprising answer will be offered by afhne iterated function systems. 

Definition 2.2. Let A be a d x d expansive integer matrix. We say that a matrix is expansive if all its 
eigenvalues have absolute value strictly greater than 1. 

Let B be a finite set of points Z d , of cardinality \B\ —: N . 

For each b G B, we define the following affine maps on M. d , 

(2.1) Tb (x) = A- 1 (x + b), {xeM d ). 

The family of functions (rb)beB is called an affine iterated function system (affine IFS). 

2.1. The attractor. Since the matrix A is expansive, there is a norm on R d such that all the maps Tb are 
contractive. Therefore the following general theorem for contractive iterated applies: 

Theorem 2.3. Hut81] Let (Ti)fL 1 be a system of N contractive maps on a complete metric space X. Then 
there is a unique compact subset K C X such that 

N 

(2.2) K=\J n (K). 

z=l 

Definition 2.4. The compact set K defined by equation (|2.2|) in Theorem 12.31 is called the attractor of the 
iterated function system (ri)f_-,. 

Example 2.5. (The middle third Cantor set) Let us consider the middle third Cantor set C. This set can 
be regarded as the attractor of an affine iterated function system in R. So d = 1. The matrix is A = 3, the 
set B is B = {0, 2}. Therefore the maps are tq(x) = § and t 2 (x) = £±2, jgR. We check the equation (|2.2p 
for the Cantor set C: To(C) is a scaled down version of the Cantor set, sitting in [0, |]; 72(C) is another 
scaled down version of the Cantor set, sitting in [|, 1]. Their union is indeed the Cantor set. 

The proof of Theorem 12.31 is based on the Banach fixed point principle: the map M : H 1— > (jfL 1 Ti(H) is 
contractive on the set of compact subsets with the Hausdorff metric. Therefore, the attractor is the limit of 
the iterations M n (H) for any compact starting set H . 

In the case of the middle third Cantor set, let us start with the unit interval [0,1]; then M([0, 1]) = 
r ([0,l])Ur 2 ([0,l]) = [0,i]U[f,l]. 

At the next step M 2 ([0, 1]) = r (M([0, 1])) U r 2 (M([0, 1])) = [0, §] U [§, \] U [§, f] U [|,1]. We see the same 
usual procedure of constructing the Cantor set. 

Remark 2.6. For a general affine iterated system (rb)beB, the equation (|2.2j) defining the attractor K =: Xb 
can be rewritten 

AX B = \J(X B + b). 

beB 

Actually, the attractor has the following representation: 

X B = ^f^A- k b k \b k G sj . 

If we assume that there is no overlap between the sets Xb + b, b G B, we can construct the function 
4>b '■— Xx B an d it satisfies 

(2.3) cf> B (A- 1 x) = J2<l>B{x-b), (x G M. d ). 

beB 
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But that is precisely a form of the scaling equation in multircsolution wavelet theory (see more below)! 
However the role of the scaling function will not be played here by the attractor of the iterated function 
system, but by an associated invariant measure. There is a direction of research which exploits equation (|2.3p . 
which enables one to construct Haar type wavelets on fractal measures. This was pursued in |DJ06a( lDMP06; 
where wavelets were constructed on Cantor sets and Sicrpinski gaskets. 

2.2. Multiresolutions. Let us recall the key facts of multiresolution wavelet theory. The purpose of mul- 
tiresolution theory is to construct wavelets, that is orthonormal bases for L 2 (M. d ) of the form 

(2.4) 1 1 det A\^ 2 ipi(A^ ■ -k) | j G Z, k e Z d , i e {1, . . . , £}} , 

where ... ,ipL are some functions in L 2 (R d ) which will be called wavelets. 
Multiresolutions were introduced by Mallat in [Mal89] . 

Definition 2.7. A multiresolution is a sequence of subspaces (V n ) n ei, of L 2 (M. d ) with the following properties: 

(i) V n C Vn+i, for all neZ; 

(ii) U„ eZ K is dense in L 2 (R d ); 

(iii) n n &V n = {0}; 

(iv) / € V n if and only if f(A-) G V n+1 ; 

(v) There exists a function (p G Vq such that 

{^(- - k)\k ez d } 

is an orthonormal basis for Vq- 
The function ip is called the scaling function. 

If one has a multiresolution, then wavelets can be constructed by considering the detail space Wq := ViQVq. 
Here on can find vectors ipi, . . . , such that their translations {ipi(- — k) \ k s Z d , i G {1, . . . , L}} form an 
orthonormal basis for Wq. Then, applying the dilation, all the detail spaces V n+ i Q V n , and the orthonormal 
basis in ( (|2.4p ) is obtained from the properties of the multiresolution. 

2.3. The scaling equation, the invariance equation. The scaling function ip in Definition ^. 7f v) satisfies 
an important equation, called the scaling equation. This is obtained by considering the function tp(A~ 1 -) 
which lies in V-\ C Vq. Since the translates of ip form a basis for Vo, the scaling equation is obtained: 

(2.5) J— < P (A- 1 x) = J2*Mx-k), (x G M d ), 

where {flk)k& is a sequence of complex numbers. 

The role of the scaling function in the harmonic analysis of fractal measures is played by an invariant 
measure. 

Theorem 2.8. Hut81 Let {ri)f =l be a contractive iterated function system on a complete metric space X. 
Letpi, . . . ,pn G [0, 1] be a list of probabilities. Then there is a unique probability measure /J, p on X such that 

N 

(2.6) Hp{E) = '^^PiUpirf 1 (E)), for all Borel subsets E. 

i=l 

Moreover, the measure /j, p is supported on the attractor of the iterated function system (Tj)£Li from Definition 

In this paper we will work only with equal probabilities Pi = . 
Definition 2.9. The unique probability measure hb satisfying 

(2.7) V>b{E) = ^(r^iE)), for all Borel subsets E, 

b£B 

is called the invariant measure associated to the affine iterated function system (Tb)beB- Equation (|2.7jl is 
called the invariance equation. 
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Example 2.10. Consider the IFS associated to the middle third Cantor set as in Example l2.5l The invariant 
measure /ib in this case is the HausdorfF measure corresponding to the Hausdorff dimension In this 
measure the piece of the Cantor set in [0, |] has measure ^, the piece of the Cantor set in [0, |] has measure 
\, etc. 

Remark 2.11. It should be stressed that the affine fractals X constructed by IFS-iterations, i.e., from the 
Cantor-Hutchinson schemes (Theorem 12.31 and Example 12. 5p . are non- linear objects. As attractors for IFSs, 
they are "chaotic"; see e.g., |Bar06j . This is similarly true for Julia sets from complex dynamics (section 1; 
i.e., invariant sets that are attractors under iteration of rational functions). They too are non- linear, and 
they carry no group structure; so there is no available Haar measure. 

We prove that Definition 12.91 offers a substitute for Haar measure in these fractal settings. One of the 
conclusions of our work (TJJOGa, DJ 07bl [JP941 IJP98| 1 is to establish that despite the chaos features of the 
attractors X, there is nonetheless a substitute for Haar measure on these sets X. Now the measures from 
Theorem 12.81 form a simplex of measures \x p indexed by p, where p = {pi)fLi ranges over the probability 
measures on {1, 2, ... , N}. Each of them plays an important role in non-commutative geometry, see |BJO04j . 

In |DJ07bj . we show that when an IFS is given, then the measures [i v on X are induced by infinite 
product measures corresponding to a choice of probability distribution (pi)fL 1 . The choice of pi — 1/N in 
Definition 12.91 is motivated by the search for spectra for our fractals X in the sense of Definition 12.11 A 
theorem in |DJ07aj states that for non-uniform weights the Hadamard condition below ( Assumption I2.2lj) 
is not satisfied, and presumably the measures hb never have spectra. (See also Problem [2]) As a result, 
the choice of measures in Definition 12.91 (i.e., uniform weights) serves as a substitute for Haar measure in a 
non-linear setting. The measures of Brolin [Bro65j supported on Julia sets in the complex plane may also 
be thought of as substitutes for Haar measure in a fractal setting where there is none. For a fixed Julia set, 
the associated Brolin measure too has the form /i p for the uniform choice of probability distribution (pt)£Lii 
i.e., the choice pi = 1/N . A Julia set X is constructed from a fixed rational function R = R(z) of degree 
N, and the idea is to get X = X(R) as an IFS attractor for the IFS resulting from a choice of N branches 
of the inverse of z i— > R{z) in the complex plane. 

We can rewrite the invariance equation for continuous compactly supported functions on M. d as follows: 

(2.8) f fdHB = ^Yl I 'f ndnB, (feC c (R d )). 

J beB J 

or 

(2.9) f f(Ax)dn B (x) = ±Yl I f{x + b)dnB{x), (feC c (R d ). 
J beB •* 



Note that both the scaling equation Ij2.5|) and the invariance equation Ij2.9|) express the dilation of an object 
in terms of the sum of translated copies of the same object. The resemblance will be more apparent when 
we take the Fourier transform of both equations. 
The Fourier transform of the scaling function is 

<p(x) = ( (fitfe' 2 ™-* dt. 



Applying the Fourier transform to the scaling equation (|2.5j) one obtains: 
Proposition 2.12. The following scaling equation is satisfied: 

(2.10) ip(A T x) — mo((A T )~ 1 x)ip((A T )~ 1 x), {xeM. d ), 
where A T is the transpose of the matrix A, and 

(2.11) meM^—L^^a^k-*, ( x(£ R d ). 

The Z d -periodic function too is called the low-pass filter in wavelet theory. 
For the invariant measure fiB, the Fourier transform is 

fi B (x) = J e 2mt d^ B {t) 

and the invariance equation implies 



s 
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Proposition 2.13. The Fourier transform of the measure hb satisfies: 

(2.12) pL B {x) ^m B {{A T y 1 x)fi B {{A T )- 1 x), (x G R d ), 
where 

(2.13) m B (x) = ^J2 e ^ ib ' X ' ( ieRd )- 

b£B 

We will call the function m B the low-pass filter of the affine IFS (r^b^B- 

We mentioned above that the wavelets are constructed from a multiresolution. The multiresolution is 
constructed from the scaling function, which in turn is built from the low-pass filter mo- If we iterate 
formally (|2.10p , and impose the condition mo(0) = 1 (from which the name "low-pass" comes from), we get 
the following: 

Proposition 2.14. The infinite product formula for if is 

oo 

(2.14) ip(x) = Y[m {{A T )- n x), (xeR d ). 

fe=i 

The infinite product is uniformly convergent on compact subsets if mo is assumed to be Lipschitz. 
In the same fashion, the Fourier transform of the invariant measure, fi B has an infinite product formula 
in terms of the low-pass filter m B (which satisfies m B (0) = 1): 

Proposition 2.15. The infinite product formula for ft B is 

oo 

(2.15) fi B (x) = l[m B ((A T )- n x), (xeR d ). 

fc=i 

The infinite product is uniformly convergent on compact subsets because is a trigonometric polyno- 
mial. 

2.4. Orthogonality of the translates of the scaling function; orthogonality of the exponentials. 

For the scaling function, one wants the translates ip(- — k), k £ J/ 1 to be orthogonal. Taking the Fourier 
transform and using a periodization, one arrives at the following condition: 

Proposition 2.16. The translates of the scaling function are orthogonal if and only if 

(2.16) h v {x) :=^\<p{x + k)\ 2 = 1 {xeR d ). 

fcez 

For the invariant measure \i B to be a spectral measure with spectrum A, we need the functions (e^)AeA 
to be orthogonal in L 2 (/j, B ). Writing the Parseval equality for the function e- x in this orthonormal basis 
one obtains: 

Proposition 2.17. The set A is a spectrum for fi B if and only if 

(2.17) h A (x) :=Y,\Mx + \)\ 2 = 1 (xeR d ). 

AeA 

Proof. The converse is also true, because, taking x = —A' for some A' <E A, the relation (|2.17p implies that 
ey is orthogonal to all the other e\. And since e- x is in the span of {e\)\^ Al the Stone- Weierstrass implies 
that this family spans the entire space L 2 (ii b ). □ 

2.5. The quadrature mirror filter condition and Hadamard triples. We saw that the scaling function 
if can be obtained from the low-pass filter mo by an infinite product formula J2.14[) . But we want the 
translates of the scaling function to be orthogonal. Combining the orthogonality condition (|2.16[) with the 
scaling equation (|2.10[) , we see that mo must satisfy the quadrature mirror filter (QMF) condition: 

(2.18) ^|m ((A T )- 1 ( 2 ; + 0| 2 -l I (ze« d ), 

where £ is a complete set of representatives for 1 d j A Tr L d . 

For the moment it is not clear what a QMF condition should be for the fractal measure fi B . One of the 
reasons is that we do not have a candidate for the spectrum A yet. However, we will see that the crucial 
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notion of Hadamard triples introduced by Jorgensen and Pedersen in [JP98j can be interpreted as a QMF 
condition. 

Definition 2.18. Let A be a dx d integer matrix. Let B, L be two finite subsets of Z d of the same cardinality 
\B\ = \L\ =: N. Then (A, B, L) is called a Hadamard triple if the matrix 



is unitary. (This is called a complex Hadamard matrix.) 

Remark 2.19. The unitarity condition in (|2.19[) may be understood as follows; d = 1. Suppose B = 
{0, bi, . . . , 6at-i}. Then unitarity in (|2 . 19[) holds for some L C K iff the complex numbers 



are roots of the polynomial 1 + z bl + z b2 + ■ ■ ■ + z hN - 1 . 

To illustrate the restriction placed on the given pair (A, B), take the example (see Remark I4.3j) when 
B = {0, 2, 3}. In that case the associated polynomial equation 1 + z 2 + z 3 = has no solutions with \z\ = 1, 
and so there is no set LcM which produces a complex Hadamard matrix. 

In the case of the middle-third-Cantor example, tq(x) = x/3 and T2(x) = (x + 2)/3; so A = 3, and 
B = {0, 2}. Up to a translation in R, the possibilities for the set L are L = {0, |(n + ^)}, for some n € Z. 
Since none of these solutions I = |(n + i) are in Z, the unitarity condition (|2.19p is not satisfied. 

Example 2.20. We will not use the example of the middle third Cantor set. It was proved in |JP98j that 
this fractal measure does not admit more than 2 mutually orthogonal exponential functions. However, in 
JP98], it was shown that a modification of this does provide an example of a fractal spectral measure. Take 
A = 4, B = {0,2}, so tqx = x/A, t^x — (x + 2)/4. The attractor of this iterated function system is the 
Cantor set obtained by dividing the unit interval in four equal pieces and keeping the first and the third 
piece, and iterating this process. 

Then one can pick L — {0, 1} to obtain the Hadamard triple (A, B, L). The unitary matrix in (|2.19[) is 



This was the first example of a (non-atomic) spectral measure which is singular with respect to the Lebesgue 
measure. 

From now on we will make the following: 

Assumption 2.21. There is set L C Z d of the same cardinality as B, \B\ = \L\ =: N such that (A,B,L) 
is a Hadamard triple. Moreover we will assume that G B and G L. 

Proposition 2.22. [LW02 Let A be a dx d expansive integer matrix, B,LcZ d with \B\ = \L\ = N. Then 
(A, B, L) is a Hadamard triple if and only if 



(2.19) 




r 2-Ki{A T y 1 (l i -l k )^ _ np 




. Based on this Jorgensen and Pedersen proved in [JP98j that the measure \ib has spectrum 




(2.20) 




where tub is the function defined in (|2.13j) . 

The proof of the proposition requires just a simple calculation. 

It is clear that we have now a QMF condition for our fractal setting. 



Definition 2.23. We say that the iterated function system 



ai(x) = {A T y 1 {x + l), (xeR d ,leL) 



is dual to the IFS (Tb)beB if (A B, L) is a Hadamard triple. 

With the Hadamard triple we have a first candidate for a spectrum of hb'- 



(2.21) 
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2.6. Lawton's theorem and transfer operators. In this section, we study a certain operator which 
is an infinite-dimensional version of a Perron- Frobcnius matrix, i.e., a matrix with positive entries. The 
Perron-Frobenius theorem applies to such matrices T, and it states that when T is given, the spectral 
radius r — r(T) is a distinguished positive eigenvalue of T. The corresponding eigenspace contains a positive 
vector. The peripheral spectrum lies on the circle centered at in the complex plane with radius r(T). There 
is an extensive theory involving a certain infinite-dimensional positive operator, which generalizes Perron- 
Frobenius in a manner which extends the Perron-Frobenius theorem from matrix theory. It is discussed in 
BJ02 , and the reader is referred to this book for details and for additional references. 

Indeed, it turns out that there is an analogous Perron-Frobenius operator in wavelet and fractal theory. We 
call it "the transfer operator" , or "the Perron- Frobenius-Ruelle operator" , but other names may reasonably 
be associated with it, see e.g., [B J02] . 

In the context of wavelets, the operator was identified in a pioneering paper |Law91| by Wayne Lawton. 
In this wavelet context, the transfer operator T is specified by a chosen filter function mo, and we denote 
the associated operator T = R mo . 

It is convenient to normalize such that the spectral radius of R mo is 1, and we will do this in our present 
paper. 

We saw that the scaling relations for wavelets (the equation (|2.5[) and for fractals (the equations (|2.2p . 
(|2.6p . and (|2.7p ) involve an identity which in a certain sense is of first order. However orthogonality relations 
involve "correlations" , and so will be of second order. But it turns out that the quantities that "measure" 
correlations, or inner product numbers for candidates of wavelet functions, or for fractal components, may 
be understood as eigenfunctions in a certain positive operator, which we call the Ruelle operator, or simply 
the transfer operator. In its present form, we give it in Definition 12. 251 below. The functions we are studying 
in Propositions 12. lBl and [2~. 171 will be cigenfunction for appropriate transfer operators. 

The best orthogonality relations in the analysis of wavelets and of fractals will come about from testing 
the leading eigenspace of the transfer operator; i.e., the eigenspace of the top eigenvalue in the spectrum of 
our transfer operator, as we make precise in Propositions 12.261 and l2~27l below. 

Suppose mo is a Lipschitz low-pass filter (m,o(0) = 1) which satisfies the QMF condition (|2.18p . Then 
one can construct the scaling function ip via the infinite product formula (|2.14p . Moreover, <p defined by 
this infinite product is a function L 2 (R d ) ( [Mal891 [Jo r06 , HW96]), and if mo is a trigonometric polynomial, 
i.e. only a finite number of the coefficients are nonzero, the infinite product will give us a compactly 
supported L 2 (R d )-function ([HW96, Theorem 3.13]). However the QMF condition is only necessary, it is 
not sufficient to guarantee the orthogonality of the translates of the scaling function. More is needed. The 
extra ingredient was provided by Lawton in |Law91| . 

Theorem 2.24. (Lawton [Law91 ) Let too be a Lipschitz low-pass filter with mo(0) = 1 and satisfying the 
QMF condition (|2.18[) . Then the scaling function (p defined by (|2.14p is an orthonormal scaling function 
(i.e., its translates are orthonormal) if and only if the only continuous 7L d -periodic functions h that satisfy 



are the constants. (Recall that C is a complete set of representatives for Z d /A T Z d .) 

Proof. The main observation needed here is that the function h v of (|2. 16|) satisfies the equation (|2.22|) . It 
can be proved continuous, since mo is Lipschitz. So if the only solutions to (|2.22p are the constants, then 
hip = 1 because h v (0) = 1. The proof of the converse statement requires a minimality property of h v : it 
is the smallest non- negative continuous solution to (|2.22p which satisfies h v {Q) = 1 (sec [BJ02, Section 5.1], 
}Dut04j ). If it is equal to 1, then there can be no other such functions. □ 

Lawton's theorem brings a new character into play: the transfer operator. 
Definition 2.25. The operator defined on Z d -periodic functions / on R d by 



(2.22) 




(2.23) 



(R mo f)(x) = K^)- 1 ^ + l))\ 2 f{{A T )-\x + I)), (x G R d ) 



is called the transfer operator. 
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In the fractal setting, one uses the function tub for the "weight" in place of mo, and the dual IFS {(Ji)i^l 
to replace the inverse branches x i— ► (A T )~ 1 (x + I). Thus the transfer operator is defined by 

(2.24) (R mB f)(x) = \m B (tTt(x))\ 2 f(at(x)), (x G R d ). 

Proposition 2.26. The following affirmations hold: 

(i) In the wavelet context, with h v defined in (|2.16|) . one has RniQ hip — hep . 

(ii) In the fractal context, if Aq is the candidate for a spectrum, defined in (|2.21[) , andh\ is the function 
defined in (12. 1T[> . one has R mB h^ — hh . 

Proof. Since we have already talked about the wavelet context (i), we only have to see what happens in 
(ii). The ideas are the same as in the wavelet context. One has to perform a computation. The important 
property needed here is that Ao = U; £ l(Z + A T Ao), and the union is disjoint (a consequence of the Hadamard 
property ([2T9]) ). For details we refer to |DJ06bj . □ 

Theorem 2.27. The set Ao in (|2.21|) is a spectrum for hb if and only if the only continuous eigenfunctions 
of the transfer operator R mB with eigenvalue 1 are the constants. 

Proof. The main ideas are the same as in the proof of Theorem 12.241 While we did not formulate this 
theorem in DJ06b|, all the required details are contained there. □ 

2.7. Cohen's orthogonality condition and the Laba-Wang theorem. A different necessary and suffi- 
cient condition for the orthogonality of the translates of the scaling function is captured in Cohen's condition. 
Before stating it, we need a definition. 

Definition 2.28. Let mo be a low-pass filter satisfying the quadrature mirror filter condition (|2.18|) . A 
finite set C := {x\, . . . ,x p -i} in [0, l) d is called a cycle if Axi = Xi+i modZ d for i G {1, ... ,p — 1} and 
Ax p = xomodZ d . The cycle C is called an mo-cycle if |mo(:Ej)| = 1 for all i G {1, . . . ,p}. 

Theorem 2.29. (Cohen Coh90a, Coh90b|) Let mo be a Lipschitz low-pass filter satisfying the QMF condi- 
tion (|2.18p . Assume in addition that mo has finitely many zeros. Let ip the scaling function defined by the 
infinite product formula (|2.14[) . Then the translates of ip are orthogonal if and only if the only mo-cycle is 
the trivial one {0}. 

Proof. The idea of the proof is to use Lawton's theorem (theorem I2.24|) and consider a continuous function 
h > such that Rm h = h. One can arrange h so that is has some zeros. If xo is a zero for h, then 
|mo((A T ) -1 (x + l))\ 2 h{{A T )- 1 {x + I)) = for all I G C. Using the QMF condition (|2T8]) . at least one 
of the terms |too((^4 t )~ 1 (x + l))\ 2 is non-zero. So one obtains another zero of h, x\ such that A T x\ — 
ccomodZ d . The process can not continue indefinitely, and the points must end in a cycle. For details see 
|Coh90al ICoh90bl ID J06bj . □ 

In the fractal context one looks for conditions when Ao defined in (|2.21[) is a spectrum for the measure 
Us- We will need a definition analogous to Definition 12.281 for iterated function systems: 

Definition 2.30. A set C := {xi, . . . ,x p } is called a cycle for the IFS (ai)i^L from Definition 12.231 if 
TiiXi = Xi+i for i G {1, . . . ,p — 1}, Ti p x p — x\ for some h, . . . , l p G L. Let ras be the function defined in 
(|2.13j) . The cycle C is called an ms-cycle if \ms{xi)\ = 1 for all i G {1, . . . ,p}. 

Theorem 2.31. (Laba-Wang LW02 ) Under Assumvtion \2.21\. suppose that ms has only finitely many zeros 
in the attractor Xl of the IFS {aiji^L- Then the set Ao defined in (|2.21[) is a spectrum for the invariant 
measure \xb if and only if the only ms-cycle is the trivial one {0}. 

Remark 2.32. Laba and Wang proved this theorem for dimension d = 1 , but their argument works also 
in higher dimensions, see [DJ06bJ. The proof follows the same ideas as in the wavelet case. 

2.8. Super-wavelets and the completion of Ao. What happens in wavelet theory if the translates of 
the scaling function are not orthogonal, so Lawton's condition in Theorem 12.241 and Cohen's condition in 
Theorem 12.291 fail? One still gets an interesting property of the associated wavelet ip in L 2 (R d ) does not 
generate an orthonormal basis but a Parseval frame in (|2.4[) (see |Law91[ lDa u92. BJQ2]). 
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Definition 2.33. A family of vectors (e,)^/ in a Hilbert space H is called a Parseval frame if 



ll/f = £l</,e<>l 2 , (/efl). 



iei 

The discovery that that the MRA wavelets are Parseval frames even in the case when they might not be 
orthogonal is usually credited to the pioneering paper of Wayne Lawton Law91 . In the book BJ02 this 
phenomenon is identified with a natural class of representations of the Cuntz algebras, one that in turn is 
motivated by signal processing in communication engineering. The authors of |BJ02j prove that the Parseval 
property follows from the axioms that defining these Cuntz algebra representations. In a more general 
context, for the generalized multiresolutions (GMRAs) the idea was extended and wor ked o ut in [BJMP05]. 

Cohen's condition implies that there are some non-trivial mo-cycles. It was shown in |Dut04| IBDP05] 
(for dimension d — 1) that these mo-cycles can be used to construct a "super- wavelet" : that is a vector ip in 
a "super-space" of the form L 2 (R) © • • • © L 2 (R) such that bay applying certain dilations and translations 
to ip (as in (|2.4p ) one obtains an orthonormal basis. Moreover, the classical wavelet ip in L 2 (R) and its 
Parseval frame is obtained by projecting the super-wavelet basis onto the first component. Also, the wavelet 
is obtained by a multiresolution technique, applied in the super-space, but with the same filter mo- 

Theorem 2.34. (Bildea-Dutkay-Picioroaga, |Dut04|, IBDP05] ) When Cohen's condition in Theorem \2.29\ is 

not satisfied, then mo generates a multiresolution orthogonal super-wavelet and an orthogonal scaling function 
for a space of the form L 2 (R) © • • • © L 2 (R). Moreover the projection of these onto the first component gives 
the classical scaling function and Parseval wavelet in L 2 (R) . 

Remark 2.35. The super- wavelet construction in L 2 (R) © • • • © i 2 (R) requires the two unitary operators: 
the dilation and translation operators. They are defined in Dut04. BDP05J. We will not include the precise 
definitions here but we will give an example which contains most of the important ideas. 

Example 2.36. (The stretched Haar wavelet) Let d = 1, A = 2 and m (x) = |(1 + e 27Tl3x ). This is a 
low-pass filter, mo(0) = 1 and it satisfies the QMF condition (|2.18|) . The infinite product formula yields 



This satisfies the scaling equation (I2.5[) . but the translates of it are clearly not orthogonal. The associated 
wavelet is 



The family in (|2.4|) is not an orthonormal basis, but it is a Parseval frame. 

Note that there is a non-trivial mo-cycle, namely {|, | }. On the Hilbert space L 2 (R) © L 2 (R) © L 2 (R) 
define the dilation and translation operators 



where U f(x) = -^f(x/2), T f(x) = f(x — 1) are the usual dilation and translation operators in L 2 (R). 
Define the super-scaling function and super-wavelet by 



1 

<P = oX[0,3)- 




U 3 (h,f 2 ,h)(x) = (Uh,Uf 3 ,Uf 2 ), r 3 (/i,/2,/ 3 ) = (Th,e 2 ^Tf 2 ,e 2 ^Tf 3 ), 



The results in |Dut041 IBDP05j prove that 

(i) The super-scaling function ip^ satisfies the scaling equation with filter mo: 

U3V3 = -j={ip3+T$ip 3 ). 

(ii) The translates of (p 3 are orthogonal: 




(iii) The super-wavelet ip 3 generates the orthonormal basis: 

{UlT^ 3 \j,k€Z} 
is an orthonormal basis for L 2 (R) © L 2 (R) © L 2 (R). 
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It is clear that results of the theory in L 2 (IR) for this example are recuperated by a projection onto the first 
component. 

Remark 2.37. The fact that the components of ips and ip^ are the same in this case is just a coincidence. 
In general, the components of the super-scaling function are defined by an infinite product formula similar 
to (|2.14p . but involving the cycle. Then the super- wavelet is defined from the super-scaling function, just as 
in the classical multiresolution theory, but using the dilation and translation of the super-space. We refer 
again to [Put 041 IBPP05] for details. 

Let us look now at the fractal measure ps, and see what happens if the condition in the Laba-Wang 
theorem is not satisfied, that is there are some non-trivial ms-cycles. In this case the exponentials e 27TlXx 
with A e A do form an orthonormal set, but the set is incomplete. It turns out that, even in this case, with 
some extra assumptions on mg, ps is a spectral measure, and each m^-cycle will contribute to a subset of 
the spectrum. Putting together these contributions one obtains the entire spectrum for ps, and Ao is just 
the contribution of the trivial m^-cycle {0}. 

Theorem 2.38. (Putkay-Jorgenscn [PJ06b ) Suppose Assumvtion \2.21\ holds andrriB has only finitely many 
zeros on the attractor of the IFS (<j;); e x,. Let A be the smallest set that contains —C for all ms-cycles C, 
and which satisfies the invariance property 

A T A + LcA. 

Then A is a spectrum for the measure p b ■ 

Example 2.39. Consider the affine IFS in Example [2T2"01 d = 1, A = 4, B = {0, 2} so t (x) = x/4, 
72(20 = (x + 2)/A. We saw that we can pick L = {0, 1} and this will yield the spectrum A = {X]fc=o | h £ 
{0,1}} for the fractal measure ps- The dual IFS is ao(x) = x/4, eri(x) = (x + l)/4. The function 
itib(x) = |(1 + e 2m2x ) will have only the trivial m^-cycle. 

But we can also pick L = {0,3}. The Assumption 12.211 is verified. The dual IFS is now ao(x) — x/4, 
o~3(x) = (x + 3)/4. There are two ms-cycles now. One is the trivial one {0} (with ao0 = 0). The other is 
{1} with CT3 1 = 1. The spectrum A3 given by Theorem 12.381 is 

As - I it AHk |fce{0,3}|u + £ I e {0, -3} 1 . 

U=0 J I fc=0 J 

Actually, one can take any L — {0,p} with p odd, and with Theorem 12.381 will yield a different spectrum 
for /j,b- And there does not seem to be any direct relation between these spectra for different values of p. 

3. Open problems 

Problem 1. f[PJ 07bj ) Under Assumption 12.211 prove that /is is a spectral measure. 

Theorem 12 .381 shows that this statement is true in dimension d = 1. In [PJ07b] we proved that this is true 
also for higher dimension if a certain extra technical condition is imposed on (A, B, L). This condition allows 
infinitely many zeros for tub in the attractor of the IFS (07)^, so the results of |PJ07bj are more general 
than Theorem 12.381 However we were not able to remove this extra condition, but we strongly believe that 
the result remains true without it. Another possibility is that the condition is automatically satisfied under 
Assumption 12.211 

Problem 2. ([LW02 ) Let /ig iP be the invariant probability measure associated to the IFS ti,(x) — p(x + f3), 
f3<E B, and probabilities fe)^, with \p\ < 1 and 8cl, \B\ = N, Pl > 0. 
Suppose ps,p is a spectral measure. Then 

(i) p = -p for some P € Z. 

(ii) Pi = ■■■ =Pn = jj- 

(iii) Suppose 6 8. Then B = aB for some a ^ and some B C Z such that there exists L that makes 
(P, B, L) a Hadamard triple. 

We obtained some results in this direction in P J07a] , but we are still far from settling this problem. 

Problem 3. Under Assumption 12.211 prove that 

^b(MXb) n 7v (x B )) = 0, (6, y eB,b^ b') 
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Thus we are interested if the afhne IFS has (measurable) overlap. 

The next Problem discusses how IFS-Hutchinson measures /i depend on the geometry in an ambient R d , 
i.e, the initialization of the IFS recursive scheme for fi. 

The related discussion of Corollary 11.11 in the Introduction involves analyticity, so the ambient space is 
the complex plane and begins with the open set, the disk {z \ \z\ < 1}. For this case, the selfsimilarity 
pattern that initializes our IFS is identified on the boundary circle. While the IFS-Hutchinson measures 
typically have a definite fractal dimension, there are many ways of arriving at them from some choice of 
transformations in an ambient M. d . 

In the next section, we discuss in detail the interplay between open sets in M. d with selfsimilarity patterns, 
and an associated class of IFS-Hutchinson measures arising from an iteration-limit applied to an initial 
geometric pattern in R d . Hence we may ask for the property in Definition ^. II (that the measure is spectral) 
to hold at the first step, i.e., for the restriction of Lebesgue measure to a suitably chosen open set ft in M. d . 
We are asking about how this property predicts the IFS-measures /z arising from a fractal iteration scheme: 
Specifically, is an IFS-measure /i spectral if its initialized counterpart is? Here "the initialized counterpart" 
refers to the data in (ii) below; and it will be taken up in further detail in the last section (Theorem 14.1 1 in 
our paper). 

We now make this theme precise in the following: 

Problem 4. Consider an affine system given by the pair (A, B) in R d with the usual conditions on the d 
by d matrix A and the subset B: 

(a) The matrix A has integral entries and is expansive; and 

(b) is in the set B, and B is contained in Z d . 
Then the following three affirmations are equivalent: 

(i) There is a subset L contained in Z d such that \L\ = \B\(= N) and (A, B, L) is a Hadamard triple. 

(ii) If the system (A, B) is initialized with an open subset f) in R d , then the restriction of Lebesgue 
measure to Q is spectral. (Is the initial measure spectral?) 

(iii) The IFS Hutchinson-measure /i constructed from the initial system (A, B) in R d by the recursive 
recipe in Definitions 12.21 and 12.91 is spectral. (Is the final measure spectral?) 

Problem 5. Characterize the polynomials z \— > R(z) for which there is an isometric isomorphism W 
of L 2 (J(R), fin) with the Hilbert space K 2 (R,P) for some finite set P(= P(R)) of distinct polynomials 
each of degree less than the degree of R, and with W taking the form F i— > (F p ) p< zp, where F(z) = 

Remark 3.1. Specifically, the requirement on the isometry L 2 (p) 3 f — > Ff S K 2 (R, P) is as follows: W is 
defined on L 2 (ii) and maps onto the subspace K 2 (R, P) in H 2 . 
Moreover 

\\f\\W) = \\Ff\\ 2 H * = 

pGP 

4. Appendix, with some of Jorgensen's recollections of conversations with Irving Segal, 

Marshall Stone, Bent Fuglede, and David Shale. 

While this paper emphasizes orthogonality, stressing fractal measures, many of the problems have a history 
beginning with the case when the measure /i under consideration is a restriction of Lebesgue measure on 
R d . The original problem was motivated by von Neumann's desire to use his Spectral Theorem on partial 
differential operators (PDOs), much like the Fourier methods had been used boundary value problems in 
ordinary differential equations. 

What follows is an oral history, tracing back some key influences in the subject, step-by-step, three or 
four mathematical generations. The reader is cautioned that it is a subjective account. And it should be 
added that analysis of fractal measures derived inspiration from a diverse variety of sources: potential theory, 
geometric measure theory, the study of singularities of solutions to partial differential equations (PDEs) , and 
dynamics. Here the focus is on operator theory and harmonic analysis, subjects which in turn derived their 
inspiration from much the same sources, and from quantum physics. A central character in our account 
below is projection valued measures as envisioned by John von Neumann vN32]and Marshall Stone Sto32]. 
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In what follows we trace influences from operator theory and mathematical physics, with projection valued 
measures playing a key role. Much of it was passed onto us from generations back. 

We begin the history of what we now call the Fuglede problem |Fug74| with conversations the second 
named author had with Irving Segal, Marshall Stone, Bent Fuglede, and David Shale. Since this never found 
its way into print, in this appendix, we draw up connections to early parts of operator theory. 

Starting in the early fifties, only a few years after the end of WWII, and still in the shadow of the 
war, mathematical physicists returned to the center of attraction around Niels Bohr in Copenhagen. Both 
Niels Bohr, and his mathematician brother Harald were still active; and still inspired members of the next 
generation. This is at the same time the generation of mathematicians who inspired me, along with a large 
number of other researchers in functional analysis, in operator theory, in partial differential equations, in 
harmonic analysis, and in mathematical physics. 

Included in the flow of scientists who visited Copenhagen around the time was some leading analysts from 
the University of Chicago, and elsewhere in the US, Marshal Stone, Irving Segal, Richard Kadison, George 
Mackey, and Arthur Wightman to mention only a few. What this group of analysts has in common is that 
the members were all much inspired by John von Neumann's view on mathematical physics and operator 
theory From conversations with Segal, I sense that von Neumann was viewed as a sort of Demi-God. 

The forties and the fifties was the Golden Age of functional analysis, the period when the spectral the- 
orem, and the spectral representation theorem, were being applied to problems in harmonic analysis, and 
in quantum field theory. Here are some examples: (1) Irving Segal's Plancherel theorem for unitary rep- 
resentations of locally compact unimodular groups |Seg50| ; (2) the Gelfand-Naimark-Segal construction 
of representations of C*-algebras from quantum mechanical states [Nai56. Nai59 ; (3) Stone's theorem for 
unitary one-parameter groups; and its generalization |Sto32| . (4) the Stone-Naimark-Ambrose-Godement 
theorem for representations of locally compact abelian groups; (5) the rigorous formulation of the Stone- 
von Neumann uniqueness theorem [vN32, vN68 ; and (6) its generalization into what became the Mackey 
machine for the study of unitary representations of semi-direct products of continuous groups, including (7) 
the imprimitivity theorem [Mac52, Mac53 for induced representations of unitary representations of locally 
compact groups. See also [SG651 iGWMl IHor55l ISha62] . 

In this context, we find much work involving unbounded operators in Hilbert space, typically non- 
commuting operators, for example axiomatic formulations of the momentum and position operators, and 
in a context where the dense domain of the operators must be taken very seriously. The simplest such 
operators are the momentum and position operators Q and P that form toy examples of quantum fields; i.e., 
the operators of Heisenberg's uncertainty inequality |GW54j . 

While the spectral theorem had already played a major role in mathematical physics, at the same time 
Lars Hormander H6r55j and Lars Garding |GW54j in Sweden were using functional analytic methods in the 
study of elliptic problems in PDE theory. 

In his first visit to Copenhagen in the early fifties, Irving Segal brought with him one of his graduate 
students, David Shale, who was supposed to have written a thesis on commuting selfadjoint extensions of 
partial differential operators, and apparently achieved some initial results; never published in any form. 

A baby example: To understand the derivative operator d/dx in a finite interval, we know that von 
Neumann's deficiency indices help. Von Neumann's deficiency indices is a pair of numbers (n, m) which serve 
as obstructions for formally hermitian operators with dense domain in Hilbert space to be selfadjoint. 

For a single formally hermitian operator in Hilbert space 7i, it is known that selfadjoint extensions exist 
if and only if n — m. The size of n (= m) is a measure of the cardinality of the set of selfadjoint extensions: 
The different selfadjoint extensions arc parameterized by partial isometries between two subspaces in Ti. each 
of dimension n. 

The property of selfadjointness is essential, as the spectral theorem doesn't apply to formally hermitian 
operators. If the usual derivative operator d/dx is viewed as an operator in L 2 (0, 1) with dense domain T> 
consisting of diffcrentiable functions vanishing at the end-points, then one checks that it is formally skew- 
hermitian, and that its deficiency indices are (1,1). As a result, we see that the derivative operator with 
zero-endpoint conditions has a one-parameter family of selfadjoint extensions; the family being indexed by 
the circle. In fact, the family has a simple geometric realization as follows: Thinking of functions in L 2 (0, 1) 
as wave functions, there is then clearly only one degree of freedom in choosing selfadjoint extensions. As we 
translate the wave function to the right, and hit the boundary point, a free phase must be assigned. This 
amounts to the 'closing up' the wave functions at the endpoints of the unit interval. What goes out at x = 1 
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must return at x — 0, possibly with a fixed phase correction. With this in mind, it was natural to ask for an 
analogue of von Neumann's deficiency indices for several operators. In fact this is a dream that was never 
realized, but one which still serves to inspire me. 

In the context of M. d , i.e., higher dimensions, it is very natural to ask the similar extension question for 
bounded open domains ft in M. d . And it was clear what to expect. 

Going to d dimensions: If functions in L 2 (ft) are translated locally in the d different coordinate directions, 
we will expect that the issue of selfadjoint extension operators should be related to the matching of phases 
on the boundary of ft, and therefore related to the tiling of M. d by translations of fl; i.e., with translations 
of ft which cover R d , and which do not overlap on sets of positive Lebesgue measure. The global motion by 
continuous translation in the d coordinate directions will be determined uniquely by the spectral theorem if 
we can find commuting selfadjoint extensions of the d partial derivative operators 




j = l,...,d, 



defined on the dense domain V of differentiable functions on ft which vanish on the boundary. We can 
take V = (ft). These d operators are commuting and formally hermitian, but not selfadjoint. In fact 
when d > 1, each of the operators has deficiency indices (oo, oo). So in each of the d coordinate directions, 
id/dxj\x> has an infinite variety of selfadjoint extensions. But experimentation with examples shows that 
'most' choices of ft will yield non-commuting selfadjoint extensions. Each operator individually does have 
selfadjoint extensions, and the question is if they can be chosen to be mutually commuting. By this we mean 
that the corresponding projection- valued spectral measures commute. 

The spectral representation for every selfadjoint extension Hj D i d/dxj\x> has the form Hj = J R A Ej (dX) , 
j = 1, . . . , d, where Ej : B (R) — * Projections (L 2 (f2)J denotes the Borel subsets of R. We say that a family 
of d selfadjoint extensions Hi , ... , H& of the respective i d/dxj operators is commuting if Ej (Aj) E k (A k ) = 
E k (A k ) E 3 (Aj) for all Aj, A k € B (R), and j^k. 

When commuting extensions exist, we form the product measure 



is an orthogonal basis in L 2 (fl). 

Remark 4.2. It might be natural to expect that an open spectral set fl will have its connected components 
spectral, or at least have features predicted by the spectrum of the bigger set. This is not so as the following 
example (due to Steen Pedersen) shows. Details below! 

If the set ft in Theorem 14. II is not assumed connected, the conclusion in Theorem 14.11 would be false. If a 
spectral set ft is disconnected, then properties of the connected components are not immediately discerned 
from knowing that there is a spectrum for fl. 

An example showing this can be constructed by taking ft = 0(p) to the following set obtained from a 
unit square, dividing it into two triangles along the main diagonal, followed by a translation of the triangle 
under the diagonal. 



E = Ei x ■ ■ ■ x E d 



on B (R d ) , and set 



where t, A S R d and A • t = J2j=i ^jtj- Then clearly 

U{t)U{t') = U{t + t'), t,t' eR d , 
i.e., U is a unitary representation of R d acting on H = L 2 (f2). 




Theorem 4.1. (Fuglede-Jorgensen-Pedersen |Fug74| [J0?8 2, Ped8?]) Let ft be a non-empty connected open 
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Details: Start with the following two open triangles making up a fixed unit-square, divided along the 
main diagonal. Now make a translation of the triangle under the diagonal by a non-zero integer amount p 
in the x-direction; i.e., by the vector (p, 0), leaving the upper triangle alone . 

Further, let 0(p) be the union of the resulting two p-separated triangles. Hence the two connected 
components in fi(p) will be two triangles; the second obtained from the first by a mirror image and a 
translation. Neither of these two disjoint open triangles is spectral; see |Fug74| . Nonetheless, as spectrum 
for fl(p) we may take the unit-lattice Z 2 . To see this, one may use a simple translation argument in L 2 (Jl(p)), 
coupled with the fact that e\(p) = 1 for all A G 1? . 

There are two interesting open connected sets in the plane that are known |Fug74| not to be spectral. 
They are the open disk and the triangle. Of those two non-spectral sets, the triangles may serve as building 
blocks for spectral sets. Not the disks! 

Remark 4.3. It is enough to assume that f£ C K d is open and has positive finite Lebesgue measure. 

If the set VI is the theorem is disconnected, then it is not necessarily true that there is a single set S, such 
that (O, S) is a spectral pair. 

To see this, take d = 1 and = (0, 1) U (2, 4), i.e., two open intervals as specified. The reason that there 
cannot be a set S such that {e s | s G S} is an ONB in L 2 (Q) is that the polynomial z A — z 2 + z — 1 does not 
have enough zeros to make (fi, S) a spectral pair. 

Proof. Fuglede proved, under more stringent assumptions, that the spectrum of E is atomic, i.e., as a mea- 
sure, E is supported on a discrete subset S C M. d . Moreover, each E ({s}) L 2 (O) for s G S is one-dimensional, 
and E ({s}) L 2 (fi) = Ce ls ' x , where e s (x) = e ls ' x is the restriction to Vt of the complex exponential func- 
tion corresponding to vector frequency s — (s\, . . . , so) G S. Since a spectral measure is orthogonal, i.e., 
E (A) E ( A') = if A R A' = 0, it follows that 

{e s \seS} 

is an orthogonal basis for L 2 (fl); in other words, we say that (f2, S) is a spectral pair. □ 

By the reasoning from the d — 1 example, by analogy, one would expect that the existence of commuting 
selfadjoint extensions will force £1 to tile R d by translations, at least if f2 is also connected. And in any case, 
one would expect that issues of spectrum and tile for bounded sets fi in M. d are related. 

Segal suggested to his U. of Chicago Ph.D. student David Shale (originally from New Zealand) that he 
should take a closer look at when an open bounded set fl in lR d has commuting selfadjoint extensions for the 
d formally hermitian, and formally commuting partial derivatives. Apparently the research group around 
Irving Segal was very active, and everyone talked to one another. At the time, Bent Fuglede was a junior 
professor in Copenhagen, and he was thinking about related matters. Moreover, he quickly showed that if 
you add an assumption, then the results for problems of selfadjoint extensions do follow closely the simple 
case of the unit interval, i.e., d = 1. 

Actually Fuglede had initially complained to Segal that the problem was hard, and Segal suggested to 
add an assumption, assuming in addition that the local action by translations in f2 is multiplicative. With 
this, Fuglede showed that for special configuration of sets fl, there are selfadjoint extensions, and that they 
are associated in a natural way with lattices L in R d . The spectrum of the representation U is a lattice L. 
By a lattice we mean a rank-ci discrete additive subgroup of M. d . 

Suppose now that d commuting selfadjoint extensions exist for some bounded open domain fl in M. d ; and 
suppose in addition the multiplicative condition is satisfied. When the spectral theorem of Stone-Naimark- 
Ambrose-Godement (the SNAG theorem) is applied to a particular choice of d associated commuting unitary 
one-parameter groups, Fuglede showed that £1 must then be a fundamental domain (also called a tile for 
translations) for the lattice L* which is dual to L. If L is a lattice in M d , the dual lattice L* is 

L* = { A G R d | A • s G 2ttZ for all s e L } . 

But more importantly, Fuglede pointed out that the several-variable variant of the spectral theorem (the 
SNAG theorem), and some potential theory, shown that if in addition £1 is assumed connected, and has 
a 'regular' boundary, then the existence of commuting selfadjoint extensions implies that L 2 (£1) has an 
orthonormal basis of complex exponentials 

{ exp(is • x) | s in some discrete subset S C 
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In a later paper J0r82 , Jorgensen suggested that a pair of sets (f2, S) be called a spectral pair, and that S 
called a spectrum of Q. With this terminology, we can state Fuglede's conjecture as follows: A measurable 
subset ft oiR d with finite positive Lebesgue measure has a spectrum, i.e., is the first part of a spectral pair, 
if and only if tiles M. d with some set of translation vectors in M. d . 

Apparently Fuglede's work was all done around 1954, and was in many ways motivated by von Neumann's 
thinking about unbounded operators. 

In the mean time, David Shale left the problem, and he ended up writing a very influential Ph.D. thesis 
on representations of the canonical anticommutation relations in the case of an infinite number of degrees of 
freedom. See Sha62j. 

Irving Segal explained to me in 1976 how his suggestions from the early fifties had been motivated by von 
Neumann. Von Neumann had hoped that all the work he and Marshal Stone had put into the axiomatic 
formulation of the spectral representation theorem would pay off on the central questions on the theory of 
PDEs. Initially, Fuglede wasn't satisfied with what he had proved, and he didn't get around to publishing 
his paper until 1974, sec |Fug74|, and only after being prompted by Segal, who had in the mean time become 
the founding editor of the Journal of Functional Analysis. 

Fuglede had apparently felt that we needed to first understand non-trivial examples that arise from tilings 
by translations with vectors that can have irregular configurations, and aren't related in any direct way to 
a lattice. Natural examples for £1 for d = 2 that come to mind are the open interior of the triangle or of the 
disk. But Fuglede proved in |Fug74| that these two planar sets do not have spectra, i.e., they do not have the 
basis property for any subset Scl 2 . Specifically, in cither case, there is no Scl 2 such that { e s |n | s G S} 
is an orthogonal basis for L 2 (ft). This of course is consistent with the spectrum-tile conjecture. 

Very importantly, in his 1974 paper, Fuglede calculated a number of instructive examples that showed 
the significance of combinatorics and of finite cyclic groups in our understanding of spectrum and tilings; 
and he made precise what is now referred to as the Fuglede conjecture f [JP94l PP00, PWMEdQ!|j0r82|). 

Fuglede's question about equivalence of the two properties (existence of orthogonal Fourier frequencies 
for a given measurable subset fl in M. d ) and the existence of a subset which makes f2 tile R d by translations, 
was for d — 1,2, and perhaps 3. But of course the question is intriguing for any value of the dimension d. In 
recent years, a number of researchers, starting with Tao, have now produced examples in higher dimensions 
giving negative answers: By increasing the dimension, it is possible to construct geometric obstructions to 
tiling which do not have spectral theoretic counterparts; and vice versa. For some of these examples we refer 
the reader to [FMM06I ITO)4l ILW95I lLW96bl ILW96cl ILW96al IlKTOSl IMatORI IKL96] . 
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